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Abstract 

We study singularity formation in spherically symmetric solutions 
of the charge-one and charge-two sector of the (2-|-l)-dimensional 5^ 
(T-model and the (4+l)-dimensional Yang-Mills model, near the adi- 
abatic limit. These equations are non-integrable, and so studies are 
performed numerically on rotationally symmetric solutions using an 
iterative finite differencing scheme that is numerically stable. We eval- 
uate the accuracy of predictions made with the geodesic approxima- 
tion. We find that the geodesic approximation is extremely accurate 
for the charge-two ci-model and the Yang-Mills model, both of which 
exhibit fast blowup. The charge-one c-model exhibits slow blowup. 
There the geodesic approximation must be modified by applying an 
infrared cutoff that depends on initial conditions. 

Mathematics Subject Classification: 35-04, 35L15, 35L70, 35Q51, 35Q60 
Physics and Astronomy Classification: 02.30. Jr, 02.60.Cb 



*This work was partially supported by the Texas Advanced Research Program, and 
carried out at the University of Texas. 

^Applied Science Fiction; 8920 Business Park Drive; Austin, TX 78759; 
jlinhart@asf.com. 

^Department of Mathematics; University of Texas at Austin; Austin, TX 78712; 
sadun@math.utexas.edu 



1 



1 Introduction, Background and Results 



The geodesic approximation is a popular tool for theorizing on the dy- 
namics of systems with known spaces of stationary solutions. As the velocity 
tends toward zero, evolution should occur close to geodesies on the space of 
the stationary solutions. 

But do equations in fact evolve as the geodesic approximation predicts? 

We study three systems to which the geodesic approximation has been 
applied. We study the charge-one sector of the (2-|-l)-dimensional 5^ cx-model 
(hereafter referred to as the "charge-one cx-model"), the charge- two sector 
of the (2+l)-dimensional S"^ a-model (the "charge-two a-model"), and the 
(4+l)-dimensional Yang-Mills model. The S'^ a-model is sometimes called 
the 0(3) (T-model or the CP^ model, and the localized solutions are called 
"lumps". In all cases the geodesic approximation suggests that localized 
solutions should shrink to zero size ( "blow up" ) in finite time. 

Blowup phenomena can be categorized as "fast" or "slow". In slow 
blowup, all relevant speeds go to zero as the singularity is approached. In fast 
blowup the relevant speeds do not go to zero. Shatah and Struwe ||T0[ proved 
that there cannot be spherically symmetric fast blowup in the charge one 
cr-model. In particular, the approach to blowup cannot be asymptotically 
self-similar, as suggested by the geodesic approximation. This theorem does 
not apply to the charge-two cr-model, or to the (4+l)-dimensional Yang-Mills 
system. 

In each case, we model the equations in a simple and robust numerical 
scheme for which the stability can be verified. Then we know the numerical 
scheme does not add features to the evolution of the differential equation. 
Consequently, we can make meaningful comparisons of the predictions of 
the geodesic approximation with the actual evolution of these equations. 
We observe slow blowup in the charge one cx-model, and fast blowup in 
the charge-two cx-model and the Yang-Mills model, in accordance with the 
Shatah-Struwe theorem. 

This is not the first paper to numerically observe blowup in these models. 
However, our intent is not merely to confirm blowup. We quantify the devi- 
ations of the true trajectories from the geodesic approximation, in particular 
for the charge-one cx-model, where these deviations are substantial. We also 
suggest a modification of the geodesic approximation for the charge-one cx- 
model, applying a dynamically generated cutoff to certain integrals to obtain 
dramatically improved accuracy. In the charge-two cx-model and the Yang- 
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Mills model, we find that the geodesic approximation is extremely accurate, 
and no such modification is needed. 

The literature on CP^ lump dynamics is extensive, especially in the 
charge-one sector, and begins with Ward's suggestion |jT^ of applying Mon- 
ton's geodesic approximation to this problem. The dynamics of lump 
scattering, in this approximation, were studied in depth by Leese ^ and 
Zakrzewski [jl5| on M^, and by Speight on S*^ |]ll| and on other Riemann 



surfaces 12 



In 1990, Leese, Peyrard and Zakrzewski studied the stabilty of numeric 
solitons, using a square grid and the traditional Runge-Kutta numerical ap- 
proach, and found that lumps on a discretized system are inherently subject 
to shrinkage. This is of direct physical interest, since many physical systems 
(e.g., an array of Heisenberg ferromagnets) are naturally discretized. How- 
ever, it raises the questions of whether a different discretization might avoid 
introducing spontaneous shrinking. In this paper we exhibit one that does. 

The results of also demonstrate the importance of doing a careful sta- 
bility analysis of whatever numerical scheme one applies, to separate proper- 
ties of the underlying differential equations from phenomena introduced by 
the discretization. The numerical methods of this paper were chosen to make 
such analysis possible [^. 

The most important numerical test of the geodesic approximation of 
blowup, for the rotationally symmetric charge-one cr-model, was done by 
Piette and Zakrzewski |^. In the charge-one sector, the Lagrangian for lump 
dynamics in the geodesic approximation has an infrared log divergence. To 
remedy this, Piette and Zakrzewski cut off their integrals at a large distance 
R, computed trajectories, and then took the R oo limit. This procedure 
predicts fast blowup (a(t) ~ to ~ ^) where a{t) is the size of the lump at time 
t), contradicting the Shatah-Struwe theorem. This prediction is also in mild 
contradiction of Piette and Zakrzewski's numerical results. They found that 
a{t) shrank slightly slower than linear, which they modeled as a power law 
with exponent slightly greater than one. 

More recently, Bizoh, Chmaj and Tabor studied the universality of blowup 
in the charge-one cr-model 0. They provide numerical evidence that the 
existence of blowup does not depend on the initial conditions (as long as there 
is enough energy), that the shape of the lump prior to blowup is universal, 
that the rate of blowup is slow {da/dt — > as t to) but that the rate is 
dependent on initial conditions. 

In this paper, we extend the results of H] in several ways. First, we obtain 
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more accurate numerical data on the rate of blowup. We see that blowup in 
the rotationally symmetric charge-one cr-model deviates from linear by a log 
correction, not a power law. Second, we show that this can be understood 
within the framework of the geodesic approximation. The cutoff Lagrangian 
of § gives the correct dynamics if we do not take the R ^ oo limit. Rather, 
R should be understood as a dynamically generated cutoff that depends on 
the initial conditions. {R is of order |a(0)/a'(0)|). The resulting trajectories 
are qualitiatively similar to those found by Speight for lump dynamics on a 
sphere |Tl|. The difference is that Speight's cutoff (the radius of the sphere) 
is independent of initial conditions, while ours is fundamentally dynamic. 

Third, we show how the shape of the lump depends on the speed of its 
shrinkage. In the geodesic approximation, the lump at all times looks like 
a static solution to the equations. We show that this is only approximately 
correct, and we quantify the deviations. In particular, we see that, shortly 
before blowup, lumps have the property that a small disk around the origin 
maps onto all of the target 5*^. 

Finally, we model the charge-two a-model and the Yang-Mills model as 
a comparison to the charge-one a-model. The form of the equations, the 
method of discretization, and the numerical stability analysis are similar for 
all three models. The conventional wisdom is that the Yang-Mills model 
should behave much like the charge-two a-model, and this is borne out by 
our numerics. Moreover, in both these cases the rate of blowup is extremely 
close to that predicted by the geodesic approximation. This gives a standard 
of accuracy by which to judge the charge-one cx-model. To wit, the deviations 
from prediction on the charge-one cr-model are much, much larger and are 
not due to numerical error. 

This paper is organized as follows. The S"^ a-model and the Yang-Mills 
model are reviewed in Section 2. In Section 3 we present our method for 
numerically integrating the PDEs, and discuss the stability and accuracy 
of this method. Our results appear in Section 4, first for the charge-one 
cr-model, then for the charge-two cr-model, and finally for the Yang-Mills 
model. 
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2 Equations being modeled 



2.1 The 2-dimensional 5^ a-Model 

The two-dimensional 5*^ cr-model has been studied extensively over the past 
few years in §, §, §,0, 0, and 



It is a good toy model for studying two-dimensional analogues of elemen- 
tary particles in the framework of classical field theory. Elementary particles 
are described by classical extended solutions of this model, called lumps. 
Such lumps can be written down explicitly in terms of rational functions. 
This model is extended to (2+1) dimensions. The previous lumps are static 
or time-independent solutions, and, with the added dimension, the dynam- 
ics of these lumps are studied. The time dependent solutions cannot be 
constructed explicitly, so in (2+l)-dimensions studies are performed numer- 
ically, and analytic estimates are made by cutting off the model outside a 
radius R. 

We will see that the 5*^ a-model displays both slow blowup and fast 
blowup. The charge-one cx-model exhibits logarithmic slow blowup, whereas 
the charge-two cr-model and the similar Yang-Mills (4+l)-dimensional model 
exhibit fast blowup. 

Identifying = CP^ = C U {oo} we can write the Lagrangian density in 
terms of a complex scalar field u: 

^„,|2 IV7„,|2 

(1) 

The calculus of variations on this Lagrangian in conjunction with inte- 
gration by parts yields the following equation of motion for the CP^ model: 



\dtu\' 




(1 + |U|2)2 


(1 + |U|2)2 



\u\ 



|2 



){d^u - dlu - dlu) = 2u{{dtuf - {d.,uf - {dyuf) (2) 



Here u represents the complex conjugate of u. 

We will look at the evolution of these equations near the space of static 



solutions. The static solutions are outlined in [Tj] among others. The entire 
space of static solutions can be broken into finite dimensional manifolds 
consisting of the harmonic maps of degree n. If n is a positive integer, then 
M„ consists of the set of all rational functions oi z = x + iy oi degree n. We 
restrict our attention to Mi and M2, the charge-one sector and the charge- 
two sector. On these two sectors the static solutions (with u{oo) finite) all 
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have the form 

M = a H (3) 

and 

respectively, depending on up to five complex parameters a, /3, 7, 5, e. To 
simplify further, we impose rotational symmetry on the solution, which is 
respected by (|^). For this regime the (2+l)-dimensional model reduces to a 
(1+1) dimensional model. These solutions become 

^ (on Ml) and ^ (on M2), (5) 

with j3 and 7 real. Note that by applying rotational symmetry we have also 
forced u{oo) to be zero.0 The geodesic approximation predicts that solutions 
evolve close to 

^ o. (6) 

z z^ 

This paper is about the actual evolution of these equations. Instead of 
looking at j3{t) or 7(t), we look for the most general rotationally symmetric 
solution to (^. We find the evolution of: 

/M) „^ /(M) p, 

Z Z^ 

The difference between the actual evolution of f{0,t) and that predicted 
in the geodesic approximation, and the difference between the profile f{r,t) 
(with t fixed) from f]{t) (which is constant with respect to r) gives us a means 
to gauge how accurately the evolution of this equation can be modeled by 
the geodesic approximation. 

It is straightforward to calculate the evolution equation for f{r,t). For 



^ There are also rotationally symmetric static solutions with u{oo) — 00. Since the equa- 
tions of motion are invariant under the transformation u — > 1/u, these behave identically 
to the solutions discussed here, and do not have to be considered separately. 
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roi u = ■> it is: 

alf . a;/ + _ + _iZ_ (,9,;). _ (Q^ff) . o. 

In both cases the static solutions are /(r, t) = constant. In the charge-one 
sector this constant is the length scale, while in the charge-two sector it is 
the square of the length scale. In the geodesic approximation, /(r, t) depends 
on t but not on r. Blowup occurs when the length scale becomes zero, that 
is f{0,t) = 0. Even apart from the geodesic approximation, a singularity 
appears when f{0,t) = 0. This is what we call the instant of blowup. 

Note that there is nothing singular about /(r, t) equaling zero when r 7^ 0. 
That merely indicates that u maps the entire circle z = re*^ to the south pole 
of S'^. We shall see that this always happens, for smaller and smaller values 
of r, shortly before blowup. 

2.2 The Yang-Mills Model 

The Yang-Mills Lagrangian in 4 dimensions is a generalization of Maxwell's 
equations in a vacuum, and is discussed at length in We can regard 
this problem as being that of a motion of a particle, where we wish our 
particles to have certain internal and external symmmetries, which give rise 
to the various geometrical objects in the problem. The states of our particles 
are given by gauge potentials or connections, denoted A, on M^, and we 
identify with the quarternions H: x = Xi + X2i + x^j + x^k. The gauge 
potentials have values in the Lie algebra of SU(2) which can be viewed as 
purely imaginary quarternions. The curvature Fij = diAj — djAi + [Ai, Aj], 
where [Ai, Aj] = AiAj — AjAi is the bracket in the Lie Algebra, gives rise to 
the potential V{A) = \Fij\'^ which is a nonlinear function of A. The static 
action is: 

S = IJ (10) 

The space of finite-action configurations breaks up into topological sectors, 
indexed by the instanton number. We are interested in the sector with instan- 
ton number one. The local minima of ( PUj) are the instantons on 4 dimensional 
space. 
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We now introduce a time variable t and consider the wave equation gen- 
erated by this potential with Lagrangian: 

L = ll (11) 

Via the calculus of variations, the evolution equation for this dynamical 
model i^ 

d^A^ = -VjFij. (12) 

Theoretical work on the validity of the geodesic approximation for an 
analogous problem, the monopole solutions to the Yang-Mills-Higgs theory 



on (3 + l)-dimensional Minkowski space, is presented in . 

The static solutions to equation (^) are simply the 4 dimensional instan- 
tons investigated in [jl|. All degree one instantons take the form: 

,, , 1 ( (x — a)dx — dx(x — a)} . . , tttt /-.on 

A{x) = - { ^ — ■ \- ^ } X = xi+X2i + x^j + Xik e e. (13) 



2\ A2 + |x - a. 
We will consider connections the form 



-4(^>^) = 2| /(^,t)+r4 r=^xl + xl + xl + xl (14) 
This is the most general form of a rotationally symmetric connection. We 



will derive an equation of motion for f{r,t) from (|12D 



To get at this, start with connections of the form 



A{r,t) = -g{r,t) {xdx — dxx} . (15) 



One can then use (|T^) to calculate the differential equation for g: 

dh = l2g'+^-^ + dlg-%gh\ (16) 



^Note that we have defined A = (Ai, yl2, Aj,. A4) to be a time-dependent connection on 
M'', not a connection on ]R^+^. However, the dynamics of Yang-Mills connections on 
is almost identical. In that model, A = {AQ^Ai) and the |9tAip term in the Lagrangian 
( inf ) is replaced by |-Fbip. After obtaining the equations of motion we apply the gauge 
choice Aq = 0. The equations of motion then reduce to (|l2|), plus a Gauss' Law constraint 
dt{diAi) = that comes from varying the Lagrangian with respect to Aq. This constraint 



is identically satisfied by the ansatz (14). Our results therefore apply to this model as well 
as to time-dependent connections on 



p4 
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Substituting g{r,t) = {f{r,t) + r^) ^, we obtain a differential equation 
for /(r,t): 

g-ff = g^j+'Al.^^ + _l_ (ojf _ (a„/)=) , (17) 

As with the charge- two sigma model, the static solutions are /(r, t) = con- 
stant, where the constant is the square of the length scale, and the geodesic 
approximation states that solutions should progress as f{r,t) = a(t). As 
before, blowup means /(0,t) = 0. 



3 Numerical Method 



A finite difference method is used to compute the evolution of these partial 
differential equations (H, |1^) numerically. Centered differences are used 
consistently except for 



d^f H and d^f H . 



(18) 



In order to avoid serious instabilities in 
special way. Let 



and 



Llf=r-\dry{dr)f = d'j + 



L2f = T~\dr)T\dr)f = dlf + 



ITD these are modeled in a 



(19) 



(20) 



These operators have negative real spectrum, hence they are stable. The 



naive central differencing scheme on ([TsD results in unbounded growth at the 
origin, but the natural differencing scheme for this operator does not. For 
the first it is 



[f{r + 5)-f{r) 



The second is done analogously. 



nV/w-/(^-^) 



(21 
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With the differencing explained, we want to derive /(r, t+At). We always 
have an initial guess at f{r,t + At). In the first time step it is f{r,t + At) = 
f{r, t) +f At with vq the initial velocity given in the problem. On subsequent 
time steps the initial guess is f{r,t + At) = 2/(r, t) — f{r,t — At). The 
appropriate forward time step can be used to compute dtf{r,t) on the right 
hand side of (||, |, |1^). Then this allows us to solve for a new /(r, t + At) 
in the differencing for the second derivative dff{r,t). This procedure can be 
iterated several times to get increasingly accurate values of f{r,t + At). In 
([T7|), for example, we iterate the equation: 



f{r,t + At) = 2f{r,t)-f{r,t-At) + {At)' 



d'rf{r,t) 



5dr f{r,t) 



2{dJ{r,t))^ 2{drf{r,t))' 8drf{r,t)r 
f{r,t) + r'^ f{r,t)+r^ f{r,t) + r'^ 



(22) 



where all derivatives on the right hand side are represented by the appropriate 
differences. 

There remains the question of boundary conditions. At the origin, in all 
models, f{r,t) is presumed to be a quadratic function, and this gives 



/(0,t) = ^/(Ar,t)-^/(2Ar,t). 



(23) 



At the r = Rmax boundary we use Neumann boundary conditions 
drf{Rmax,t) = 0. We essentially have a free wave equation outside of this 
boundary in all of the equations. This boundary condition allows informa- 
tion that has propagated out to it to bounce back toward the origin, but it 
does not permit information to come in through the boundary from infinity. 

In the models generated from (^, results indicate that the appropriate 
form for /(r, t) at any given time t, is a parabola in r, rather than a straight 
line. When looking at this aspect of the model, the f{Rmax,t) boundary 
condition was changed to respect the parabola: 



drfi^Rmaxjt) dj-f(^Rmax Ar,t) 



Rr, 



R„ 



Ar' 



(24) 



3.1 Stability and Accuracy 

A lengthy analysis of the stability of the numerical schemes for (|], |1^) is 
performed in 0. That for (|^) should be directly analogous to that for (p^. 
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We do not wish to go into the full details of the analysis for the stability or 
convergence in this paper, as they are tedious. 

However, we do emphasize the importance of such an analysis. 

Discrete numerical models, such as these, are usually approximations to 
a desired continuum model. Usually the true continuum solution is only 
obtained as some parameter goes to zero or to infinity. In evaluating a model 
numerically, we introduce error. We have the unavoidable round off error 
that comes with every computer model, and also a truncation error that 
occurs from only taking finitely many terms. 

Sometimes a method that seems otherwise favorable has undesired be- 
havior that is caused by the discretization used in the problem. 

For example, questions arise from |Q about the stability of the lumps for 
(H). They found that a soliton would shrink without perturbation from what 
should be the resting state. This seems to be a function of the numerical 
scheme used for those experiments. The numerical scheme used here for 
equations (||, 0) has no such stability problems. Experiments show that 
a stationary solution is indeed stationary unless perturbed by the addition 
of an initial velocity. 

In unstable numerical schemes, the roundoff error can come to be com- 
mingled with the calculation at an early stage and can then be magnified 
until it overcomes the true answer. A good, basic discussion of stability and 
numerical schemes can be found in ||^. 

Apart from the stability analysis, one can measure the effect of the spa- 
tial discretization by tracking different trajectories with identical initial con- 
ditions, a small fixed At, and a varying Ar. Any scheme will break down 
when the system develops important features that are too small compared to 
Ar. When two trajectories start to differ, that implies that the computation 
with the larger value of Ar has started to break down. 

We see this in Figure |l], where we track f{0,t) for the charge-one S'^ a- 
model. We pick input parameters of /(r, 0) = fo = 1.0, dtf{r,0) = vo = 
-0.01, and At = 0.001. We track evolution with Ar = 0.1, Ar = 0.05 
and Ar = 0.01. The difference between f{0,t)/\r=o.i and /(O, t)AT-=o.oi 
at the point of blowup (specifically, when /(O, t)Ar=o.oi hits zero) is less 
than 0.01. This is considerably less than the grid size difference of 0.09. 
Furthermore, the difference when /(0,t)Ar=o.i = 0.1 is smaller still, with 
1/(0, t)Ar=o.i — /(O, t)Ar=o.oi| ~ 0.002. This correspondence in the three tra- 
jectories indicates that the influence of the grid size on the features of the 
model can be ignored until /(0,t) has shrunk roughly to size Ar. 
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f(0,t) with delt = 0.001 , deir varying, closeup of blowup 



0.25 



0.15 



0.1 



0.05 















del 
deir 
deir 


= 0.1 

= 0.01 

= 0.05 































































85 90 95 100 105 110 115 

Time 

Figure 1: Different trajectories for different Ar values in the charge-one a- 
model. We have zoomed in on a neighborhood of the blowup time in order 
to be able to see where these trajectories begin to deviate from one another. 

4 Comparison to the Geodesic Approxima- 
tion 

We wish to compare the actual evolution of the blowup of the charge-one 
and charge-two S"^ cr-models and the Yang-Mills model with that predicted 
by the geodesic approximation. 

In our ansatz for all three models, we model the evolution of /(r, i), and in 
all cases the geodesies on the moduli space of stationary solutions are /(r, t) = 
a{t). When we calculate deviations from the geodesic approximation of this 
profile, we calculate deviations between the actual shape /(r, T), with T 
fixed, and f{r,T) ~ a. Note also that outside of a sufficiently large ball 
about the origin, all three equations are well approximated by the linear 
wave equation dff = df.f . The interesting behavior should occur near the 
origin, and so evolution toward blowup is tracked as the evolution of /(O, t). 
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4.1 Charge-One 5^ a-Model 



For the charge-one S"^ a- mo del, our calculations follow those done in M and 
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In IP the time evolution of blowup via the shrinking of lumps was studied. 
The Lagrangian was cut off outside of a ball of radius R, to prevent logarith- 
mic divergence of the integral for the kinetic energy, and then an analysis of 
what happens in the R oo limit was performed. In |jll| the problem of 
the logarithmic divergence in the kinetic energy integral is solved by inves- 
tigating the model on the sphere S"^. The radius of the sphere determines a 
parameter for the size analogous to the parameter R for the size of the ball 
the Lagrangian is evaluated on in 

Equation gives us the Lagrangian for the general version of this prob- 
lem. In the spherically symmetric geodesic approximation for the charge-one 
sector, we consider functions of the form 



u 



z 



(25) 



If we restrict the Lagrangian to this space we get an effective Lagrangian. 
The integral of the spatial derivatives of u gives a constant, the Bogomol'nyi 
bound, and hence can be ignored. If one integrates the kinetic term over the 
entire plane, one sees it diverges logarithmically, so if /3 is a function of time, 
the solution has infinite energy. 

Nonetheless, this is what we wish to investigate. We presume that the 
evolution takes place in a ball around the origin of size R. Up to a multi- 
plicative constant, the effective or cutoff Lagrangian becomes 



' \r^ + P? 



which integrates to 



{dtff 



In 1 + -TT 



For geodesies, energy is conserved, so 



{dtff 



In 1 + 



P 



P 



R' 



P + R\ 



R^ 
P + R^ 



c 
~2 



(26) 



(27) 



(28) 
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f(0,t) vs ct with the best fit line to the data 
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Figure 2: Predicted plot of /(O, t) vs ct for the cutoff Lagrangian with R 
100 in the charge-one a-model. 



with c (and hence c^/2) a constant. Solving for dtf we obtain 

c 



dtf 



In 1 + 



R' 



R' 



PJ P + R' 



(29) 



Since we are starting at some value /o and evolving toward the singularity 
at / = this gives: 



cdt. 



(30) 



In the R oo limit, as taken in 0, this gives linear evolution, i.e., fast 
blowup. But this would violate the Shatah-Struwe theorem [|l^]. However, 
taking a finite value of R gives a/| ln(t)| corrections and slow blowup, in 
accordance with the theorem. 
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Evolution of Origin 




.0.1 I 1 1 1 1 1 1 1 

50 60 70 80 90 100 110 120 

Time 

Figure 3: Plot of f{0,t) vs t for an actual trajectory in the charge-one a- 
model. The slight nonlinearity is reminiscent of the prediction shown in 
Figure |[ 

In equation (|30|) , the integral on the right gives ct. The integral on the 
left can be evaluated numerically for given values of R, fo and /(0,t). A 
plot can then be generated for ct vs. /(0,t). What we really are concerned 
with is f{0,t) vs. t, but once the value of c is determined this can be easily 
obtained. One such plot with /q = 1.0, R = 100 of f{0,t) vs ct is given in 
Figure |^. This curve is not quite linear, as seen by comparison with the best 
fit line; both are plotted in Figure |^. 

The scale invariance of this problem allows us to consistently take /o = 1.0 
without loss of generality, and the key dimensionless parameter is dtf{0,0). 
We compare this with the numerical computer model, which was run with 
various small velocities. The initial velocity is dtf{r, 0) = Vq. Another impor- 
tant input parameter is Rmax, the largest value of r modeled, which is chosen 
large enough so that information from the origin cannot propagate out to the 
boundary and bounce back to the origin within the time to blowup. Other 
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1/[fdot(0,t)]'^2 vs ln[f(0,t)] and best fit line to the data 
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Figure 4: Actual plot of l/{dtf{0,t)Y vs ln(/(0,t)) for the charge-one a- 
model and the best fit line to this data. The cutoff Lagrangian predicts 
linear dependence, while the R ^ oo limit predicts that I /dtp should be 
constant. 



input parameters are Ar and At. 

We track /(0,t) as it heads toward this singularity, and find that its 
trajectory is not quite linear, as seen in Figure ^. This is suggestive of the 
result obtained in the predictions for this model. 

We want to check the legitimacy of the prediction ([30| ) , pictured in Figure 
^ This requires a determination of the cutofiQ R and the kinetic energy 
We already have /o and /(O, t). To determine R and c, observe from equation 

■^The Lagrangian cutoff R is not to be confused witli Rmax, the maximum value of r 
modeled numerically, i? is a parameter of the geodesic approximation, not of the numerical 
integration. 
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f(0,t) vs. t: Model and Predicted 




Figure 5: Comparison of the cutoff geodesic approximation to actual data in 
the charge-one a-model. The predicted course of f{0,t) vs. t from equation 
(|30| ) is practically indistinguishable from the actual trajectory. 



that: 



1 



In 1 + 



Since R is large and / is small, we can rewrite equation (^TJ) as 

1 MR') - ln(/^) - 1] 



(31) 



(32) 



The plot of ln(/) = ln(/(0,t)) vs l/idtff = l/(9t/(0, t))^ should be linear 
with the slope = and the intercept = (21n(i?) — l)/c^. Such a plot 
is easily obtained from the data, and the parameters c and R are readily 
calculated. 

Figure || is a plot of ln(/(0,t)) vs. 1/(9^/(0, t))^, with initial conditions 
Ar = 0.01, At = 0.001, /o = 1.0 and Vq = -0.01. If the correct dynamics 
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Table 1: Best-fit parameters c and i? as a function of vq, with /o = 1.0 for 
the charge-one cr-modeL 



^^0 


c 


R 


-0.005 


0.0145 


115 


-0.00667 


0.0187 


89 


-0.01 


0.0263 


53 


-0.0133 


0.0342 


49 


-0.02 


0.0485 


34 


-0.03 


0.0683 


25 


-0.05 


0.104 


17 


-0.06 


0.121 


15 



(33) 



were given by the R oo hmit of equation (^), this graph would be of a 
horizontal line, since the velocity would be unchanging. On the other hand, 
if equation (jS^) is correct with R finite, we should get a line with nonzero 
slope. It is easily seen that the plot of ln(/(0,t)) vs. l/{dtf{0,t)Y is nearly 
a linear relationship, as predicted by the cutoff Lagrangian, but is not quite 
a straight line. This may indicate that the effective values of R and c are 
themselves changing slowly with time. 

The best fit line to this data has slope -2810 and intercept 10200, corre- 
sponding to c = 0.0267 and R = 62.1. Using these values of c and R in the 
calculation of equation (|30|), we obtain the plot of /(O, t) vs t given in Figure 
^ This is overlaid with the model data for /(0,t) vs. t for comparison. 
These two are virtually indistinguishable. 

We have seen that modifying the geodesic approximation by applying a 
fixed and finite cutoff greatly improves its accuracy. Moreover, corrections to 
linear behavior are seen to be logarithmic, not power-law. However, we have 
not yet demonstrated any predictive power, since the cutoff was deduced 
from the data only after the fact. What is needed is an understanding of 
how the dynamically generated cutoff R (and the kinetic energy c^/2) vary 
with initial conditions. 

Table |l| gives the best-fit values for c and i? as a function of the initial 
velocity Vq, under the initial conditions /o = 1.0, Ar = 0.01 and At = 0.001. 
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R vs 1/vO and best linear fit to the data 



R vs 1/vO data 
y = 0.5407*x + 6.032 



20 40 60 



100 120 140 160 1£ 

1/vO 



200 



Figure 6: The dynamical cutoff R grows linearly with 1/vq in the charge-one 
(j-model. 



R is roughly linear in 1/fo, with best-fit line 

^ + 6.032. (34) 

This fit is shown in Figure A linear fit with 1/vq makes sense, since the 
time to blowup is itself of order I/vq. Ris roughly half the radius of the light 
cone at t = 0. 

Next we consider the shape of the lump as we approach blowup. That 
is, what is /(r, T) with T fixed? The geodesic approximation indicates that 
f{r,t) should be independent of r, but this is only approximately true. The 
graph of f{r,t) versus r stays close to horizontal, although there is some 
slope downward as time increases. This is shown in Figure |^. 

Making a closer inspection of the time profiles /(r, T), as in Figure |^, we 
observe that the initial part of the data is close to a hyperbola as seen in 
Figure |[ This fit is close, but not nearly as close as those for profiles in the 
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Time Profiles of Data 
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20 
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100 



Figure 7: Time slices /(r, T) for the charge-one a-model. 



charge-two sector of the S"^ cr-model or the (4+l)-dimensional Yang-Mills 
Model, below. While the asymptotic conditions for the hyperbola are not 
consistent with those of a solution to this equation, this form need only be 
assumed on a ball about the origin, outside of which the function is of the 
form f{r>R,T)= a. 



4.2 The Charge-Two cr-Model 

Equation (Q) gives us the Lagrangian for the general version of this problem 
We are using 

X{t) 



u 



(35) 



for our evolution, and via the geodesic approximation we restrict the La- 
grangian integral to this space, to give an effective Lagrangian, as we did 
in the previous section. The integral of the spatial derivatives of u gives a 
constant, and hence can be ignored. Unlike the charge-one cr-model, this 
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Time Profile f(r,72) and best hyperboiic fit to tfie data 

0.33 I 1 1 1 




0.27 - 

0.26 - \- 

0.25 I 1 1 1 1 

20 40 60 80 100 

r 

Figure 8: Time slices for the charge-one cr-model evolve a hyperbolic bump 
at the origin. The transition from a hyperbola to a horizontal line occurs at 



Lagrangian does not have an infrared divergence, and there is no need to 
apply a cutoff. Up to a multiplicative constant, the effective Lagrangian is 

By conservation of energy, 

dtX = 2ciVA, (37) 
for some constant ci, and we integrate to obtain 

\{t) = {c,t + C2)^ (38) 

which we rewrite as 

X{t)^p{t-tof. (39) 
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Table 2: Parameters for best fit parabola to f{0,t) vs. initial data /o and Vq 
for the charge-two S*^ a- mo del. 



/o ^^0 to p 

1.0 -0.01 200 0.0000250 

1.0 -0.02 100 0.0000998 

1.0 -0.03 67 0.000224 

1.0 -0.04 50 0.000398 (^0) 

1.0 ^ 100^2 0.0000499 

1.0 ^ 200^2 0.0000125 

1.0 ^ 200^3 0.00000833 



Equation ( pUD is how we predict that f{0,t) will evolve. 

As before, the scale invariance of the model allows us to fix /o, and the 
other input parameters are the same as before. A typical evolution of /(O, t) 
looks like that given in Figure ^. (That figure is actually a trajectory for 
the Yang-Mills model, but the results from these two models are essentially 
the same.) In every case, the trajectory of /(O, t) fits a parabola of the form 
pit — toY almost perfectly. The dependence of the parabolic parameters 
p and to on initial velocity Vq is given in Table 0. In all cases we took 
/(r,0) = /o = 1.0, dtfir,0) = vq, Ar = 0.025 and At = 0.001, and in all 
cases we obtain 

2 v'^ 
to = — and p= (41) 

l^^ol 4 

The predictability of the parameters from Table ^, and the closeness of the 
parabolic fits indicate there is no correction needed to the predictions of the 
geodesic approximation. 

With the evolution of /(O, t) taken care of, we consider the shape of the 
time slices f{r,T) for a given fixed T. The geodesic approximation suggests 
that the graph of f{r,T), which starts horizontal, should remain horizontal. 
Instead, an elliptical bump forms at the origin, as seen in Figure 0. (That 
figure is also from the Yang-Mills model, which exhibits essentially identical 
behavior) . 
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If we describe the elliptical bump by the equation 



- + - - 1 (42) 

then the parameters a, b and k evolve as 

a = t (43) 

b = (44) 

k = 1.0 + vot. (45) 

The growth of this ellipse suggests that the curve is trying to adjust from a 
horizontal line to a parabola 

f{r,t)=pr^ + h, (46) 

where 

To test this, we began runs with initial data 

/(r, 0) = 1 - vy/8; dtf{r, 0) = -Vq. (48) 

for various values of vq. The time slices of these runs retained a parabolic 

profile, as shown (for the Yang-Mills model) in Figure |TT|. That is, we expect 



/(r,t)--|r^ + |(t 



2 

1^0 1 



(49) 



Substituting this closed- form expression into (|]), we find the error is 

±,.'-'^r^(t-^)\ (60) 

64 32 V \vo\J 

Within a fixed radius, and for small fo, this error remains small for all time. 
However, the parabolic fit cannot be extended out arbitrarily far, especially 
insofar as the profile ( ^6]) has logarithmically divergent potential energy. 
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Evolution of the Origin 




-0.05 I 1 1 1 1 1 1 1 1 

40 60 80 100 120 140 160 180 200 

Tiitie 

Figure 9: Sample trajectory of /(O, t) vs. t in the Yang-Mills model. The 
actual trajectory is almost exactly parabolic. 

4.3 The (4+l)-Dimensional Yang-Mills Model 

Equation (pJ]) gives us the Lagrangian for the general version of this problem. 
We are using 

^=2\ f + r^ ] r= ^xl + xl + xl + xl (51) 

and we apply the geodesic approximation. Restricting the Lagrangian to the 
moduli space gives us an effective Lagrangian. The portion of the integral 
given by 

(52) 

represents the potential energy and integrates to a topological constant, hence 
it may be ignored, as in the study of the 5^ a-model. We need to calculate 

L = \l \dtA,?d'x. = [ %^d'x = const, x (53) 
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Time Profiles of Data 
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Figure 10: Time slices /(r, T) for the Yang-Mills model evolve an elliptical 
bump at the origin. The edge of the bump is always at r = T. 



Since energy is conserved, dtf is proportional to f^^'^, and we integrate to 
get 

f=p{t-tof. (54) 

In short, the predicted evolution is identical to that of the charge two cr- 
model. 

A typical evolution of /(0,t) is given in Figure ^. In this figure, the 
graph of 0.000025()f: — 200)^ neatly overlays the graph of /(O, t). This picture 
represents the evolution where /q = 1.0 and Vq = —0.01. 

If we now look at the evolution of the initial line, /(r, 0) = /o, we obtain 
the same striking result as with the charge-two 5*^ cx-model. The initial line, 
/(r, 0) = /o, evolves an elliptical bump at the origin that grows as time passes. 
Figure |10| shows this behavior. The clear indication here is that corrections 
to the profile are quadratic. The elliptical bumps can be modeled exactly 
as before. Also as before, this suggests that we look for an approximate 
solution of the form f{r,t) = p{t)r'^ + h{t). Calculating from the general 
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Time profiles witin parabolic initial data 
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Figure 11: With parabolic initial conditions, the time slices in the (4+1)- 
dimensional Yang-Mills model do not form bumps, but instead remain 
parabolic. 



form of our ellipse in (^2|) we get the same results for p as in p7|), and h 
should evolve the same way as /(O, t). When a run is started with this initial 

data dth = vq = -0.01, /(r, 0) = 1.0 - pr^ and p = -f = -0.0000125, the 



time slices of the data retain a parabolic profile. This is shown in Figure |Tl 
The curvature of the parabola at the origin, as measured by the parameter 
p from equation ( PB| ) changes by less than 1 part in 100 during the course of 
this evolution. 

The general form of a parabolic /(r, t) is the same as in equation (^Uj). If 
we substitute this into the partial differential equation ([T7|), we get an error 
of 

2 



r\ (55) 

Since our concern is the adiabatic limit, fo is always chosen to be small, and 
this difference goes to zero quickly in the adiabatic limit. 
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The form of the connection A given by this parabohc model is 




Mr.t) = -{- — -^). (56) 



This is close to the connection A given by the geodesic approximation mul- 
tiplied by an overall factor. This form suggests a relativistic correction. 



5 Conclusions 

We have studied the time evolution of solutions towards a singularity in the 
charge-one 5*^ cx-model, the charge-two 5^ a-model and the (4+l)-dimensional 
Yang- Mills model. We use a simple and robust numerical scheme to do so, for 
which the stability has been verified. The importance of the stability analy- 
sis is that we have no problem trying to separate real effects from numerical 
artifacts. 

In the charge-one a-model, the geodesic approximation is ill-defined due 



to the divergence of the Lagrangian. A theorem of |ITO| tells us that there is 
no fast blowup in this model. The divergence of the geodesic approximation 
is solved by cutting off the Lagrangian at a radius R < oo. In the R ^ oo 
limit of this cut-off Lagrangian, we have a prediction for linear evolution, 
which would be fast blowup in violation of the theorem in [T^. However, 



the cut-off Lagrangian itself is extremely accurate in predicting the shrink- 
ing of lumps toward the singularity in this model. This evolution is close 
to linear with complicated logarithmic corrections (equation (|30|) ) going ap- 
proximately as A/ln(to — t). Corrections from linear are not power-law, nor 
power-law corrected by a logarithm, as suggested in [§]. 

We model the charge-two cr-model and the (4+l)-dimensional Yang-Mills 
model to provide a comparison to our charge-one S"^ a-model results. We 
verify the belief that the charge two S"^ cr-model is similar to the (4+1)- 
dimensional Yang-Mills model; the numerical results show that the two mod- 
els are almost identical in behavior. There is no theorem restricting fast 
blowup in either of these models, and we find fast blowup as predicted by 
the geodesic approximation in both. The numerical models' evolution toward 
the singularity is in lock step with prediction, and no corrections are needed 
to this order. Since the stability analysis of the numerical schemes in the 
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charge-one cr-model and the Yang-MiUs model is similar, the accuracy 
of the numerics to prediction in the charge-two S'^ cr-model and Yang-Mills 
model provide a accuracy scale on which to judge the charge-one S"^ cr-model. 

In addition to modeling the evolution toward the singularity, we also can 
characterize the shape of the profiles of the lumps with the time variable fixed. 
In the geodesic approximation, with the rotationally symmetric functions we 
have chosen to model, the predicted evolution in all three cases was /(r, t) — 
a(t). The actual profiles in all three cases have quadratic corrections. 
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